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Abstract. In this paper, we consider the Frenkel-Kontorova model of a one 
dimensional chain of atoms submitted to a potential. This potential splits 
into an interaction potential and a potential induced by an underlying sub- 
strate which is a quasicrystal. Under standard hypotheses, we show that every 
minimal configuration has a rotation number, that the rotation number varies 
continuously with the minimal configuration, and that every non negative real 
number is the rotation number of a minimal configuration. This generalizes 
well known results obtained by S. Aubry and P.Y. le Daeron in the case of a 
crystalline substrate. 
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1. Introduction 

The Frenkel-Kontorova model |FK| describes the physical situation of a layer of 
a material over a substrate of other material (see for instance |BKp . In the one 
dimensional case, the layer of material is described by the configurations of a bi- 
infinite chain of particles on the real line. These configurations are parametrized by 
a bi-infinite non decreasing sequence (^n)„gz ^^^^ numbers, where 0„ represents 
the position of the particle labeled by n. 
The potential energy of the chain reads: 

where U describes the interaction between particles (only interactions with the 
nearest neighbors are considered), and is a potential induced by the substrate 
and depends on its nature. 

The following standard extra asumptions are made on U and V: 

• Smoothness: the functions U and : M — » R are C^; 

• Convexity: U"{x) > 0, Va; G M; 

U(x) 

• Behavior at oo: lim — — — — +oo. 

Even if the above sum is only formal, it is possible to look for equilibrium con- 
figurations which minimize locally the energy (ground states). More precisely let 
us consider the function 7i : M x R ^ M defined by: 

n{e,9') = u{e-9') + v{e). 

For a configuration (^n)„g2, let us set: 

j=p-i 
j=o 

We say that the segment {9i, ^i+i, . . . , Oi+p) of the configuration {On)nez minimal 
if 

'Hp{9i,9i+i, . . . ,6i+p) < TLp{9[,6[^i, . . . ,6[j^p), 

for any other segment (0-, . . . , 6'-_,_p) such that 9[ = 9i and 6'-^p = A 
configuration {9n)n£z minimal if all its segments are minimal. 

The substrate is a crystal when the configuration of the chain of atoms it is made 
of, is an increasing sequence QC = {sn)nez, such that there exists q G and L > 
verifying: 

Sn+q = s„ + L, V n e Z. 

In this case it is natural to consider that a potential V associated with the crystal 
QC is a periodic C"'-function with period L: 

V{9 + L) = V{9), V 61 e R. 

This situation when the substrate potential is periodic has been described by S. 
Aubry and P. Y. Le Dearon. Their seminal work |AD| . together with the indepen- 
dent approach of J. Mather gave rise to the so called Auhry- Mather theory, 
which yields in particular a good understanding of minimal configurations. 
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Let /5 £ M, a configuration (6'n)„gz ^ rotation number equal to p if the limit: 

lim — = P- 

n^±oo 71 

Let us remark that the inverse of the rotation number can be interpreted as a 
particle density. 

Aubry and le Daeron proved in particular that any minimal configuration has a 
rotation number, that the rotation number is a continuous function when defined 
on the set of minimal configurations equipped with the product topology, and that 
any positive real number is the rotation number for some minimal configuration^. 

The aim of this paper is to consider the case when the substrate is a quasicrystal 
in order to derive, in this more general setting, a similar description of the set of 
minimal configurations. 

To fix notations and definitions, let us consider a bi-infinite substrate chain of 
atoms represented by its configuration (sn)„gx- Two segments (s„, . . . , s„+p) and 
{sq, . . . , Sg+p) are said equivalent if there exists r G M such that: 

Sq+i = Sn+i + T, Vi = 0, . . . ,p. 

The chain QC = (sn)^^^ is a quasicrystal if the following properties are satisfied^(see 
for instance ILPp : 

• Finite local complexity 

For any M > 0, the chain possesses only finitely many equivalence classes 
of segments with diameters smaller than M. 

• Repetitivity 

For any segment S in the chain, there exists R > such that any ball with 
radius R contains a segment equivalent to S. 

• Uniform pattern distribution 

For any segment S in the chain, and for any point a; S M, the quantity 

n{S, X, M) 
M 

converges when M — > +cxd uniformly in x to a limit that does not 

depend on x, where n{S, x, M) denotes the number of segments equivalent 
to S in the interval [x,x + M]. 

Notice that a crystal (with period L) is a quasicrystal and in this particular case, 

for each segment S in QC, one has: 

^iS) = ^, 

where p(S') stands for the number of segments equivalent to S* in a period L. 

For any i? > 0, a function Vqq : M ^ M is a potential with range R associated 
with a quasicrystal QC if for each pair of points x and y in R such that 

QCDBnix) - x^ QCnBR{y) - y, 

we have: 

Vqc{x) = VQciy), 



^Actually Aubry-Mather theory says much more about the combinatorics of minimal configu- 
rations when projected on a circle with length L. 

^See Proposition 12. l| for a dynamical interpretation. 
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Figure 1 . Construction of the Fibonacci chain 

where Bm{z) stands for the baU with center z and radius M. Whenever QC is a 
crystal with period L, it is clear that a potential with range R > associated with 
this crystal is a periodic potential with period L. 

We call short range potential associated with a quasicrystal QC a potential with 
range R, for some R > 0. 

Example: A standard example of quasicrystal is given by the Fibonacci sequence. 
Consider the set G of configurations (s„)^ such that: 

• sq is located at 0; 

• the lengths of the intervals [sn,s„_|_i] have two possible sizes: either large 
and equal to L or small and equal to S. 

The substitution: 

( L ^ LS 
\ S ^ L 

induces a map '5 on t/ defined as follows: 

For a sequence (s„)„ in Q, consider the sequence of lengths G {L, S}°° defined 
by In = Sn+i — Sn, '^ri G Z. Applying to each In the substitution rule we get a new 
sequence (0„ G {L, S}°° . The new configuration {s'n)n — ^((■Sn)„) is obtained by 
setting: 

• s'o = 0; ^ 

• S'n+i = s'n+l'n, Vn G Z. 

Starting with the equidistributed configuration (sJ'J^, where Sn+i — Sn — L, Vn G Z, 
it is easy to check that the sequence of configuration converges when 

k — * +0O (for the product topology) to a configuration (s^)„ . This configuration 
is on the one hand a quasicrystal and on the other hand a periodic point with period 
2 of the operator 5*. This quasicrystal is called the Fibonacci chain (See Figure^. 

There are several ways to construct a short range potential associated with 
the Fibonacci chain. A simple one consists in choosing two real valued smooth 
functions, and vs,l with compact support on the interval {^1,1) where 

< 2/ < S{< L). A potential Vpib with range 2L, can be defined as follows 
(see Figure 121): 

• for each n G Z and for each G (sjf - /, + /): 
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Figure 2. A short range potential associated with the Fibonacci chain 

- Vpibid) = vl^l{Q-s'^) if both intervals [sjf.i, s^] and [sjf , sj^^] have 
the same length L; 

- VFib{0) = vsx{9 - sjf ) if the intervals [s,'f„i,s^] and [sjf ,s^i] have 
different lengths. 

• iovOi U„ez(.<. - I. s'n + I). VF^b{e) - 0. 
The main result of this paper is the following theorem: 

Theorem 1.1. ^ 

For the Frenkel-Kontorova model with a short range potential associated with a 
quasicrystal: 

(i) any minimal configuration has a rotation number; 

(ii) the rotation number is a continuous Junction when defined on the set of 
minimal configurations equipped with the product topology; 

(iii) for any p > 0, there exists a minimal configuration with rotation number p. 

It turns out that, once the appropriate objects have been defined, the proof of 
Theorem 11.11 has the same structure as the modern proof for crystals that can be 
found for instance in jB] or More precisely, in the crystal periodic poten- 

tial factorizes through a real valued function defined on a circle. In the quasicrystal 
case, a short range potential factorizes through a real valued function defined on a 
more sophisticated compact metric space called the hull of the quasicrystal. This 
hull possesses locally the product structure of an interval by a Cantor set i.e it is 
a solenoid. This solenoid can be seen as the suspension of the action of a minimal 
homeomorphism on the Cantor set. 

Minimal homeomorphisms on the Cantor set have been extensively studied in 
topological dynamics and possess a powerful combinatorial description in terms 
of Kakutani-Rohlin towers (see for instance |CPSj ). The aim of Section |21 which 
is devoted to the substrate, is to rephrase these well known results in our specific 
context, namely for a suspension, in order to see the hull as an inverse limit of one 
dimensional branched manifolds. These branched manifolds will play a central role 
in the proof. 

In the crystal case, when projecting a minimal configuration on the circle, the 
Aubry-Mather theory shows that it wraps around the circle in a very special way, 
namely it is ordered as the orbit of a degree one homeomorphism of circle. In 
the quasicrystal case, there exists also some combinatorial obstructions, they are 
described and analyzed in Section 13 which is devoted to the ground states of the 
overlying layer. 

■^Prom a more physical point of view, it is straightforward but interesting to rephrase Theorem 
ll.ll in terms of particle density of minimal configurations. 
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Section 0] is devoted to the proof of Theorem 11.11 First, as for the crystal case, 
we show, using the inverse hmit structure of the hull given in Section El and the 
combinatorial obstructions gotten in Section |31 that minimal configurations have 
a rotation number (point (i)). Then we prove (again as in the crystal case) the 
continuity of the rotation number (point (ii)). The proof of point (Hi) of Theorem 
ll.ll in the crystal case is done first by constructing periodic minimal configurations 
for any positive rational rotation number and then to use the continuity of the 
rotation number to get a minimal configuration for any prescribed positive rotation 
number. In the quasicrystal case, the scheme is exactly the same, but the set of 
rational numbers needs to be replaced by another dense subset of the positive reals. 
More precisely when the rotation number is not 0, its inverse has to be a finite 
linear combination with positive integer coefficients of the densities of patches of 
the quasicrystal. 

This paper ends with two final remarks developed in Section |S1 the first one con- 
cerns dynamical systems. In the case of a crystal, minimal configurations for the 
Frenkel-Kontorova model are orbits of a twist map on an open annulus. Similarly, 
in the quasicrystal case, these minimal configurations are also orbits of a dynamical 
system that we describe. The second one consists in giving the bases of a possible 
extension of the theory to quasicrystals in higher dimension. 

Remark: It should be pointed out that one can find in the literature several studies 
on the the Frcnkel-Kontorova model with a quasi-periodic potential, for instance a 
potential which is the sum of two periodic potentials with incommensurable periods 
(see for instance ' EFR.Jp . Such potential cannot arise naturally from an underly- 
ing one dimensional substrate. Actually, the underlying object which organizes the 
minimal configurations and which was a circle in the crystal case and a solenoid in 
the quasicrystal case, becomes a 2-torus. More precisely the real line is immersed as 
a line with irrational slope in the 2-torus. Actually, this is a situation more complex 
than the one we are dealing with in this paper which essentially uses dimension 1 
objects, and this explains the lack of exact results in this quasi-periodic case. 

Nota Bene: To avoid an unnecessary dichotomy and unless explicitly specified, the 
quasicrystals that appear in the sequel will not be crystals. 

2. The hull of a quasicrystal 

In this section, we recall some background results concerning quasicrystals. Most 
of these results are true in any dimension and they are presented here in the par- 
ticular case of the dimension 1. Material for Subsections 12 . II and 12 . 21 can be found 
in |KPj . |BBG' and |BG| . For Subsections 12.41 and 12.51 a discrete approach can be 
found in [HPS] and we refer again to (BBG) . ,BG) and (Sj for a more geometrical 
point of view. 

2.1. The hull as a dynamical system. Consider a chain of atoms whose con- 
figuration is a quasicrystal QC = (sn)„gz- clear that each translated copy 
QC — u = {sn — u G I^, of QC is again a quasicrystal. 

The set of translated copies QC+R of a quasicrystal can be equipped with a topology 
that, roughly speaking, says that two quasicrystal configurations are close one to 
the other if in a big ball centered at in R, the segments of both configurations 
inside the ball are equivalent and equal up to a small translation. Such a topology 
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is metrizable and an associated metric can be defined as follows (see |RW| for more 
details): 

Consider two quasicrystal configurations QC — u\ and QC — U2 in QC + R. Let A 
denote the set of e e]0, 1[ for which there exists u with |m| < e, such that QC — ui 
and QC ~ U2 + u coincide in Bi/^{0). Then 

S{QC - ui, QC - U2) = inf A if ^ ^ 

<5(QC-ui, QC-U2) ^ 1 if yl = . 

Hence the diameter of QC + M is bounded by 1 and the M-action on QC + M is 
continuous. The continuous hull Q{QC) of the quasicrystal QC is the completion of 
the metric space {QC + M.,S). 

As a direct consequence of the finite local complexity property, it is easy to check 
(see for instance |RWp that 17 (QC) is a compact metric space and that any element 
in Q{QC) is a quasicrystal whose segments are equivalent to segments in QC. The 
translation group M acts on fl{QC) and the dynamical system {Q{QC), M) possesses 
(by construction) a dense orbit (namely the orbit QC + R). On the one hand, the 
repetitivity property is equivalent to the minimality of the action i.e all its orbits 
are dense, (see |KPp and, on the other hand, the uniform pattern distribution is 
equivalent to the unique ergodicity i.e the M-action possesses a unique invariant 
probability measure (see BG ). These results yield the following proposition. 

Proposition 2.1. Let QC be a quasicrystal, then the dynamical system (il(QC),R) 
is minimal and uniquely ergodic. 

In the sequel, we will denote by fi the unique probability measure on fl(QC) which 
is invariant under the K-action. 

2.2. The canonical transversal. The canonical transversal, rio(QC), of the hull 
^{QC) of a quasicrystal QC is the collection of quasicrystals in il{QC) which contain 
{i.e. such that one atom in the chain is located at 0). 

Proposition 2.2. (see |KP| ) The canonical transversal of a quasicrystal is either 
a finite set when QC is a crystal or a Cantor set when not. 

It follows that when the quasicrystal QC is a crystal, n{QC) is homeomorphic to a 
circle and when not Q{QC) has a solenoidal structure, i.e. it is locally the product 
of a Cantor set by an interval. 

The return time function C : no{QC) is defined by: 

C{T) = inf{< >0\T - t e noiQC)} VT e ^^{QC). 

The finite local complexity implies that the function C is locally constant, it takes 
finitely many distinct values Li, . . . , Lp and the clopen (closed open) sets Ci — 
C~^{Li) for i = 1, . . . ,p form a partition of floiQCY (see FigureOJ. 
The first return map r : Qq{QC) fio(QC) is defined by: 

t{T) = t - c{T) vr e n{QC). 

The unique invariant probability measure /i of the R-action on f2(QC) induces a 
finite measure v on Qq{QC) which is r- invariant. 



Recall that clopen sets form a countable basis for the topology of a totally disconnected set. 
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Figure 3. The time return function 

For any i = 1, . . . ,p and for any clopen set C in Cj, the measure v satisfies: 

HC) = ^K{{T-u) Tec, ue [0,LJ}). 

The subsets of il{QC) which read C — u where C is a clopen set in one of the Ci's 
and u G [0, Li[ are called verticals. 

The following lemma is a direct byproduct of the above definition: 

Lemma 2.3. For any S > 0, there exists a positive constant eQc{S) such that, for 
any vertical V with diameter smaller that eQciS) and any pair of configurations 
QC — X and QC — y in V, we have: 

QC n Bs{x) - x = QCn Bs{y) - y. 

2.3. Potentials on the hull. The following result shows that a short range po- 
tential associated with a quasicrystal QC, factorizes through a function on ri(QC). 

Lemma 2.4. Let QC be a quasicrystal, and let Vqc be a continuous short range 
potential associated with QC. Then, there exists a unique continuous function Vqc : 
f^(QC) R such that: 

VQcix) = VQciQC-x), Vx e R. 

Furthermore, when Vqc has range R > 0, there exists a positive constant eQc{R) 
.such that Vqc is constant on each vertical with diameter .smaller than eQc{R). 

Remark: Notice that when QC is a crystal, Lemma |^] simply means that for any 
continuous periodic function g : M — > R with period L, there exists a continuous 
function G : M./L.'L R such that g = Gon, where tt : R R/i.Z is the standard 
projection. 

Proof of Lemma \2.4\ Assume that Vqc is a potential with range i? > 0. Applying 
Lemma 12.31 for any vertical V with diameter smaller than eQc{R) and any pair 
QC — x and QC — y in V , we have: 

QC n Br{x) - x^ QC n BR{y) - y, 

and thus: 

Vqc{x) = VQciy). 
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Since the set QC + M fl y is dense in V, it follows that a continuous function Vqc 
which satisfies Vgc (2;) — Vqc{QC ~ x), Va; G M, must be constant on and equal 
to VQc{y) for any real number y such that QC — y € V. Conversely the function Vqc 
defined this way is clearly continuous, satisfies Vqc{x) — Vqc{QC — x), W x £ K, 
and is constant on verticals with diameters smaller than eQc{R)- O 

2.4. Kakutani-Rohlin towers. The following construction, which has been de- 
veloped for the study of minimal dynamics on the Cantor set, will be useful all 
along this paper. It is often referred to as Kakutani-Rohlin towers (see |HPS| ). 
Choose S > and fix a clopen set C in one of the Ci 's with diameter smaller than 
eQciS). 

Consider the first return time function Cc associated with this clopen set (which 
is constructed exactly as the first return time function in rio(QC)). The finite 
local complexity hypothesis implies that the function Cc is locally constant and 
takes finitely many values Lc^i, ■ ■ ■ ,Lc^p(c)- The clopen sets I?c,i = 'C^^(^c,i) for 
i = 1, . . . ,p{C) form a partition of C. Again because of the finite local complexity 
hypothesis, there exists a finite partition of C in clopen sets £j, j = I, . . .r such that 
for each j S {1, . . . , r}, there exists i G {1, . . -piC)} so that the following properties 
are satisfied: 

• £j C Dc,i; 

• for each u G [0, Lc^i[, £j — u is a vertical with diameter smaller that eQc{S). 
For j = 1, . . . , r, the set: 

{Sj - u, G [O.LcA}, 

is called a tower with height Lc^i- The union of all these towers realizes a partition 
of f2(QC) and the data (QC, S,C, is called a Kakutani-Rohlin towers 

system with size S. 

For j = 1, . . . , r, consider the set £j C Dc^i and for each u G [0, Lc.i[, we call 
floor of the tower £j x [0, ic,i[i the vertical £j — u. By identifying all the points in 
this vertical, each tower projects on a semi-open interval and the whole hull Q.{QC) 
projects onto a smooth branched one-dimensional manifold which is a collection of 
r of circles 71, . . . , 7^ tangent at a single point. This branched manifold is called the 
skeleton of the Kakutani-Rohlin tower system (QC, S,C, {£j}ji^{i,,,,_r})- It inherits 
a natural orientation, a differentiable structure and a natural metric respectively 
issued from the orientation, the differentiable structure and the Euclidean metric 
of the real line R (see Figure 0J. We denote it B and call tt : Vl{QC) B the above 
identification. 

The proof of the following lemma is plain. 

Lemma 2.5. Let R > and Vqc a continuous potential associated with QC with 
range R > 0. Consider a Kakutani-Rohlin towers system with size S and let B be 
its skeleton. 

Assume that S > R, then the function Vqc : r2(QC) —>■ M induced by Vqc descends 
to a continuous function Vqc : B 

Vqc on ^ Vqc ■ 

Whenever the function Vqc is -smooth for some < r < 00, then the function 
Vqc is also -smooth. 
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Figure 4. A towers system and its skeleton 



2.5. Inverse limits. Let us choose an increasing sequence (S'„)„>q going to +00 
with n and let us construct inductively an infinite sequence of Kakutani-Rohlin 
towers system as follows (see |HPSp : 

• Fix a point xq in flolQC). 

• Choose a clopen set Co containing xq, with diameter smaller than egc(»S'o) 
and construct a Kakutani-Rohlin towers system (QC, So, Co, {^Oj"}je{i,...,ro}) 
with size 5*0 . Up to a renaming of the indices, we can assume that xq 
belongs to £^0,1 ■ We denote by Bq the corresponding skeleton and call 
TTo : 0(QC) — > Bq the standard projection. 

• We choose a clopen set Ci C £0,1 which contains xq with a diameter 
small enough so that we can construct a Kakutani-Rohlin towers system 
(QC, Si,Ci, ri}) with size S'l such that each of its towers inter- 
sects all the towers of the previous system. Up to a renaming of the indices, 
we can assume that xq belongs to £1.1. We denote by Bi the corresponding 
skeleton and call vri : il(QC) Bi the standard projection. 

• Assume we have constructed a sequence of nested clopen sets C„ C C„_i C 
. . . Ci C Co containing xq and, for each p — 0, . . . , n, a Kakutani-Rohlin 
towers system {QC, Sp^Cp, {£p,j}j^[i ^.p}) with size Sp such that each of 
its towers intersects all the towers of the system associated with p — 1, and 
such that xo belongs to Sp^i . We iterate the procedure by choosing a clopen 
set Cn+i C £n,i which contains xq small enough so that we can construct a 
Kakutani-Rohlin towers system (QC, S'„+i, C„+i, {£'n+i j}je{i,...,r„+i}) with 
size Sn+i such that each of its towers intersects all the towers of the system 
associated with n. Up to a renaming of the indices, we can assume that xq 
belongs to £n+i.i- We denote by Bn+i the corresponding skeleton and call 
7r„ : fl{QC) Bn the standard projection. 

For each n > 0, fix a point y in Bn+i- The set T^n+iiu) included in a floor of a 

tower of the tower system (QC, Sp, Cp, {Spj}je{i r }), and thus descends through 

7r„ to a single point on Bn- We have defined this way a continuous surjection: 

Tn ■ Bn+1 Bn- 
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The inverse limit: 

lim Bn = {{xn)n>o I ^" ^ ^'^'^ r„(x„+i) = x„, Vn > 0}, 
gives a re-interpretation of the hull fl{QC): 

Proposition 2.6. (BG., When equipped with the product topology the set lim Bn is 

homeomorphic to n{QC). 

Notice that the map r„ : S„+i Bn induces a p„ x Pn+i homology matrix 
Mn whose integer coefficient rn„.ij is the number of times the loop Jn+i,j in Bn+i 
covers the loop jn,i of Bn under the action of the map t„. We remark that the 
construction of the sequences of towers systems we made insures that, for all n > 0, 
the matrix Af„ has positive coefficients. These matrices carry information about 
the invariant measure v on the Cantor set through the following lemma (see for 
instance [UEHI): 

Lemma 2.7. 

J^nA = X! Vi G {1, . . . ,_p(n)}, 

where Vn.i is the measure of the clopen set £n,i- 

Again the following lemma is plain: 

Lemma 2.8. Let R > and Vqq be a continuous potential associated with QC 
with range R > and choose an increasing sequence {Sn)n>a Qoing to +oo, such 
that R < Sq. Then, for each n > 0, the function Vqc induces on each branched 
manifold Bn a function VQCn which satisfies: 

yQC,n °Tn = Vgcri + l- 

3. Combinatorics of minimal configurations 

In this section, we consider the minimal segments for a short range potential 
with range R associated with QC. 

Lemma 3.1. Let L and J = I + u be two disjoint intervals in M such that for each 
e in I: 

Br{9) nQC + u = Br{6 + m) n QC, 
and let {6i, . . . , On) be a minimal segment such that [9i, On] contains I and J . For 
any pair of consecutive atoms 0^ andO^+i inlClQC, the interval [Om + u,Om+i+u] 
contains at most two atoms of the minimal segment. 

Proof. The proof works by contradiction. Assume that there exists a pair of atoms 
Om and Om+i in / H QC, such that the interval [Om + u, O^+i + u] contains three 
consecutive atoms of the minimal segment, say Oi, 0/+i, and Oi^2- 

[0l,0l+2] C [0,n+U,0„,+i+u]. 

We consider the new segment obtained by taking the atom in position and 
assigning to it the new position Oi+i — u (Figure [SJ. When u > (what we can 
assume without loss of generality) this segment reads: 

{Oi, . . . ,0i, . . . , Om, Ol+i — U, Om+l, . . . , 0l,0l+2, ■ ■ ■ , On)- 
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To get a contradiction wc arc going to show that the potential energy of this new 
segment is smahcr than the potential energy of the first one. On the one hand, 
since -Bh(0;+i) n QC — u = Bfi{9i+i — u) Ci QC, the potential energy induced by 
the substrate on the atom that changed its position, keeps the same value: 

VQci0i+i) = Vqc{9i+i - u). 

Thus, the sum of the potential energy induced by the substrate on the whole seg- 
ment is not affected by this change of position. 

On the other hand, the difference of the potential energy of interaction between the 
new segment and the former one is given by: 

Au = {uiOm - 01+1 + u) + u{ei+i -u- 6m+i) - u{em - em+i)) 
-{u{Oi - ei+i) + u{ei+i - 61+2) - u{ei - 01+2)). 

Let us introduce the new variables: 

X = 6m- Oi+i +u, Y = Oi+i -u- Om+l, 
X' = 61-61+1, ¥' = 61+1-61+2. 



We have: 
and: 



X < X' <Q and Y < Y' < Q, 



AU = {U{X) + U{Y) - U{X + Y)) - {U{X') + U{Y') - U{X' + Y')). 
For t e [0, 1], let us consider the function: 

G{t) = U{tX + {l-t)X') + U{tY+{l-t)Y') -U{t{X + Y) + {l-t){X' + Y')). 
We have: 

AU = G(l) - G(0), 

and 

G'{t) = U'{tX + (1 - t)X'){X - X') + U'{tY + (1 - t)Y'){Y - Y') 
- U'{t{X + y) + (1 - t){X' + Y')){X + Y-X'-Y') 
= {U'{tX + (1 - t)X') - U'{t{X + Y) + {1- t){X' + Y'))){X ~ X') 
+ {U'{tY + (1 - t)Y') - U'{t{X + y) + (1 - t){X' + Y'))){Y - Y'). 
Observe that for t e[0,l]: 

tX + {l- t)X' > t{X + r) + (1 - t){X' + Y') 

and 

tY + {l- t)Y' > t{X + Y) + {1- t){X' + Y'). 

Using the convexity of U, more precisely the fact that U' is an increasing function 
we get that: 

AU < 0, 

and this inequality is strict as long as 6m ^ 6i — u and 6m+i 7^ 61+2 — u. In this 
case, wc get the desired contradiction. 

In the situation when 6m = 61 — u and dm+i = 61+2 — u, we remark that both 
segments {6m, 61+1 — u,6 — m+1) and {61-1,61,61+2) are not minimal and thus 
the new configuration wc constructed is not minimal. The corresponding minimal 
segment (by fixing the extremities 6*1 and On) has an energy which is strictly smaller, 
a contradiction. □ 
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Figure 5. Move of a single atom in a segment 

The following lemma shows that there are actually more obstructions than the 
ones described in Lemma l3.il 

Lemma 3.2. With the same hypotheses and notations as in Lemma \3.1\ consider 
two disjoint pairs of successive atoms 9m < ^m+i < ^m' < Om'+i in I d QC, such 
that at least one of the four points 9m + u < 9m+i + u < 9m' + u < 9m'+i + u does 
not belong to the minimal segment. Concerning the two intervals [9m + u^9m+i-'ru\ 
and [9m' 9m'+i +u\, none of the following three situations is possible (see Figure 
0); 

(i) both intervals contain two atoms of the minimal segment; 

(ii) both intervals do not contain atoms of the minimal segment in their interi- 
ors; 

(iii) one of the interval contains two atoms of the minimal segment and the other 
does not contain atoms in its interior. 

Proof. As for Lemma 13.11 we are going to reach a contradiction assuming that 
situation (i) occurs. The proof for the other two cases works exactly along the 
same lines. Let 9i < 9i+i < 9i' < 6';'+i be atoms of the minimal segment such that: 

[9i,9i+i] c [9 

and 

[6l' , 6l' + l] C [9m' + U, 9m' + 1 + u] . 

Assuming again that u > 0, let us move some atoms of the minimal configuration 
to reach the following new configuration: 

{9i, . . . , 9m,9i + i — U,. . . ,9l' — U, 9m' + l, ■ ■ ■ , 9i,9m+l + U, . . . 9m' + U, 9l' + i, . . . , 9n)- 

Since for each 9 m I: 

Br{9) r\QC + u = Br{9 + u)r\ QC, 

the potential energy induced by the substrate on the atoms did not change even 
if the atoms have changed their positions. Thus, the sum of the potential energy 
induced by the substrate on the whole segment is not affected by this change of 
position. 

On the other hand, the difference of the potential energy of interaction between the 
new segment and the old one is given by: 

MJ = AUi + AC/2, 

where 

AUi = iUi9m - 9i + i +u) + U{9i - 9m+i - u)) - {U{9m - 9m+l) + U{9i - 9l+i)), 
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Figure 6. The forbidden 3 situations 



and 



AC/2 = {u{ei, -u-em'+i)+u{em'+u-ei,+i))-{u{em' -em'+i)+uiei, -ei^+i)). 

Let us introduce the new variables: 

Xo = 9m Xi = 9i+i -u and Yq = 9i - u Yi = 9m+i- 

We have: 

AUi = ([/(Xo - Xi) + U{Yo - Yi)) - ([/(Xo - Fi) + U{Yo - X^)). 
This yields: 

AUi = - / U"{v - u)du I dv. 

Jxo \Jxi 



Since U is convex, Xq < Yq and Xi < Yi and at least one of these inequalities is 
strict, we get: 

AUi < 0, 

and for the same reason 

AU2 < 0. 

This yields a contradiction. □ 

From the previous two lemmas, we deduce that the quantity of atoms of the 
minimal segments which belong to / and to I + u differ by an integer smaller than 
2. This is summarized in the following proposition that will be our main tool in 

the sequel of this paper. 

Proposition 3.3. Let {9i, . . . ,9n) be a minimal segment and let I be an interval 
in [9i,9n], then there exists an integer N £ Z'^ such that for any pair of disjoint 
intervals Ii = I + u\ and I2 = I + U2 in [9i , 9n] which satisfy that for each 9 in I 
and k = 1, 2: 

Br{9) nQC + Uk = Br{9 + Uk) n QC, 
each interval Ik contains either N, N + 1 or N + 2 atovfis of the minimal segment. 
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4. Proof of Theorem 11.11 

4.1. Existence of a rotation number. In this subsection, we consider a mini- 
mal configuration for a potential with range R associated with QC. Let us con- 
sider an increasing sequence (Si)^^^ going to -l-oo with I and such that Sq > 
R and consider also an associated sequence of Kakutani-Rohlin towers systems 
(QC, Si,Ci, j }je{i....,r,});>Q and the corresponding sequence of skeletons ('B„)„>q 
as constructed in Subsection 12.51 
The identification 

I : X eM.^ QC ~ X E n{QC) 
induces an immersion of the real line in QC and the image of a configuration 
through this immersion is an element (^n)„ in Q{QC)^ where 0„ = QC — 9n, for 
all ri G Z . In turn, for any I > 0, the projection ni : fl{X) — » Bi transforms this 
sequence in an element in Bf where = ni{9n) — t^i o I{9n), for all n E Z. 

Furthermore we have: 

VQcA^n) = Vsc(^rO = VQc{0n). 

The following lemma is a direct consequence of Proposition 13. 31 
Lemma 4.1. Let be a minimal configuration such that lim 0„ = -t-oo and 

lim On = — oo (resp. let {6p, . . . , 6q) be a minimal segment). Then, for any I > 

n — ^ — oo 

and any j G {1, . . . ,n}; there exists an integer Nij such that for each loopjij ofBi, 
each connected component of {tti o I)~^{'yi,j) C M (resp. each connected component 
of (tt; °I)~^{'Ji,j) C M which does not intersect (— oo, 6p] U [6q, +oo) ) contains either 
Nij or Ni_j + I or Nij + 2 atoms of the minimal configuration (resp. the minimal 
segment). 

In other words, when n increases, the projection of the minimal configuration (resp. 
the minimal segment) stays the same amount of time in a given loop up to an error 
of 2. 

Now we can prove the existence of a non negative rotation number for any 
minimal configurations. 

First, consider a minimal configuration (0„) such that lim 0„ = -|-oo and 

lim 9n — — oo. Let us estimate the length of the interval [6*0, ^n] for n > 0. Let 

n — ^ — oo 

nij be the number of times tt; o /([Sg, 6'„]) covers completely the loop -fij of Bi. We 
have, for each I > 0: 

pi PI 
'^'"■i.jLi^j < On-Oo < '^n^jLij + 2Li, 

where Lij is the height of the tower associated with the loop jij {i.e. the length 
of the loop ji j) and 

Li = max Li , . 
3e{i,...,pi} 

On the other hand we have: 



pi pi 

^«;jiV,,, < n < J2ni,j{NLj+2) + 2(7Vj + 2), 
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where 



This yields: 



Ni = max Ni ,-. 
je{i,...,pi} 



pi pi 



Pi — „ — Pi 

E n,,,(7V,,,- + 2) + 2(Ar; + 2) ^ ^i.i^^^J 

j=i i=i 

When n goes to +oo the quantity: 



pi 



goes to the measure of the clopen set Ei^j. It follows that the sequence {dn — 
9o)/n has bounded limit sup and limit inf and that any accumulation point p of 
this sequence satisfies: 

Pi Pi 

E '^ij^ij E '^ijLij 

< p < 



PI — — Pi 

Recall that the measure u is the transverse measure associated with an invariant 
probability measure on the hull ri(QC) and thus: 

pi 

^"ijLij = 1. 



On the other hand we have: 



Pi 



J2 ^^'i = 

We deduce that: 

1 1 

< p < . 

Pi — 1^ — pi 

3=1 i=i 
Since these last inequalities are true for any ^ > 0, and since ^{Ci) goes to as Z 
goes to +00, it follows that the sequence {On — Oo)/n converges to the limit: 

lim W- 

E ^i.j^ij 
i=i 

Observe that this rotation number is difi^erent from 0. 

Consider now a minimal configuration which satisfies lim On = M < +oo or 

n—*+oo 

lim On = m < +00. The constant configuration 

n—*+oo 

On = Oo, Vn e Z, 

has obviously a rotation number equal to 0. Let us assume now that the minimal 
configuration is not constant and satisfies lim On = M < +oo. Let us show that 
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we cannot have hm 9n = —oo. Indeed, consider the interval [M — 2R, M + 2R] 

n — > — oo 

and choose u > such that the interval [M — 2R — u,M + 2R — u] is disjoint from 
[M - 2R, M + 2R] and such that: 

B2r{m -u)nQC + u = B2r{M) n QC. 

Consider now, for n large enough, the interval [6'„ — -R, 0„] C [M — 2R,M]. The 
number of atoms in [0„ — i?, 0„] goes to +00 with n. If lim 0„ —00, it follows 

n — ^ — 00 

from Proposition 13. 31 that the number of atoms in [0„ — i? — u, On — u] and thus in 
\M — 2R — u,AI + 2R — u] , goes to +co with n. Consequently the minimal sequence 
(^n)ji has an accumulation point in [M — 2R — u,M + 2R — u] when n goes to —00 
which is a contradiction. Thus for a minimal configuration we have: 



lim 6n < +00 

n — >-\-oo 



lim 6n > —00 

n — > — 00 

{On)n is bounded 



{On)n has rotation number 0. 
This ends the proof of Part {i) of Theorem ll.il 

4.2. Continuity of the rotation number. Consider a sequence {Om^n)^ of niin- 
imal configurations with rotation numbers p„i which converges, in the product 
topology, to a minimal configuration (6'„)„ with rotation number p > 0. We fix 
I > and choose a loop 7;^ in Bi. Consider the first time when, starting from on 
the real line and going in the positive direction, the projection of the configuration 
{9n)n enters in this loop. Let us do the same for the configuration (0m,n)„- Since 
{(^m,n)n converge to (6'n)„ in the product topology, for m large enough both config- 
urations stay the same time in the loop for their first visits. It follows from Lemma 
14. II that the minimal number of times Nm.i.j, the projections of the configurations 
{6m,n)n Spend in the loop "/ij of Bi, and the minimal number of times Nij., the 
projection of the configuration {0m,n)n spends in the same loop 7;^ , satisfy: 

\Nm,i,j - Ni^,\ < 2, yj e {!,..., pil)}. 

The rotation number p,„ of the configuration {0m,n)„ satisfies: 

1 



< 



< 



1 



2iyiCi) 



On the other hand 



1 



pi 



< o < 



J2 Vl,j^m,l,j 



1 



2v{Ci) 



pi 



This implies that for m large enough: 

1 1 

P Pm 

Considering bigger and bigger I yields: 



< 



HCi)- 



lim Pm = p. 

m — >+oo 
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When the rotation number p = 0, we have proved that the configuration {On)„ is 
bounded. Let M be its upper bound and consider the loop 70.2 in Bq on which M 
descends by projection. If M falls on the singular point, we consider the loop where 
the M — e's for e > small enough, are falling. Fix K > 2, when m is big enough, 
the projection of the configuration (0m,n)„ (whose rotation number is assumed to 
be different from 0) must spend at least a time K in the loop jo,i during one of its 
visits and thus, thanks to Lemma f4. II at least K — 2 times at each of its visits. It 
follows that the rotation number of {Om.n)^ satisfies: 

1 

(K - 2)1^0^^ 

and, consequently: 

lim = 0. 

m — ^+00 

Thus, we have proved Part (m) of Theorem ll.il 

4.3. Construction of minimal configurations. Observe that a constant config- 
uration is a minimal configuration with rotation number 0. For positive rotation 
numbers, we are first going to construct minimal configurations for a dense subset 
of rotation numbers in M+. 

The good candidate to be a dense set in for which minimal configurations 
can be construct is suggested by the expression (*) obtained in the previous sub- 
section. Again, let us consider an increasing sequence (5'i)i>o goi'^S to +cxd with I 
and such that 5*0 > R- Consider also an associated sequence of Kakutani-Rohlin 
towers systems (QC,S';,C;,{f;j}jg|i_ ,,!});>o ^^^'^ corresponding sequence of 
skeletons {Bn)n>o constructed in Subsection l2.5l RecalHng that the i^ij's are the 
measures of the clopen sets £ij, we define the set J- as follows: 



^-<(i5r^ , VTVi,, e z+ \ {0}, Vj e {i,...,pz}, V? > o 

Since the measures of the clopen sets £ij go to zero with I uniformly in j, we check 
easily that !F is a dense subset of K"*" . 

Proposition 4.2. For any real number po in T , there exists a minimal configura- 
tion with rotation number po- 

Proof. Fix Iq > and choose pi^ positive integers Ni^^i, . . . , Ni^^^p^^ . Consider the 
positive real number: 

= -0 ^ 

J=l 

Let us construct a minimal configuration with rotation number po- 
Step 1: For j = l,...,p/g, consider on the loop 7;„j of the oriented branched 
manifold Bi„, Ni^j — 1 points 6/0,1 < • • • < big,Ni j-i, disjoint from the singular 
point ^[^{xo) of Big (where we recall that n;>o C; ~ {^a}) as shown in Figure [7| 
For each j — 1, . . . ,pig, we consider the segment: 



(7rio(a;o), 6/0,1, . . . ,5/o,Ar,^,^-i,7r/o(a;o)). 
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Figure 7. The branched manifold with its marked points 

Thanks to Lemma 12.81 it makes sense to compute the potential energy of this 
segment and to consider the position of the points, fo;o,i, ■ • ■ , big^Ni^ j-i, which min- 
imizes this potential energy. Having done it for all loops, we denote Bi^ the col- 
lection of these marked points (all the fo/o,fc's and 7r/Q(a;o)) on Big and consider the 
subset of the real line {nig o I)~^{Big). It is a discrete subset that we can ordered 
as a bi-infinite increasing sequence (6'/o,ti)„- This subset of R contains the subset 
{nig o I)~^{TTig{xo)) which is a quasicrystal. The configuration {Oig^n)^ is made with 
a concatenation of minimal segments whose extremities are consecutive points in 
(tt/q o /)^^(7r;„ (xq)) and, there are exactly pig different equivalence classes of seg- 
ments, each of them corresponding to a minimal segment starting at the beginning 
and ending at the end of a loop in Big . 

Step 2: Consider now the subset Ti^^{Big) of the branched manifold Big+i. This 
subset contains the singular point 7r;Q+i(a::o) and for each j = the 
loop 7/q+i J of Big-i-i contains Nig^ij — 1 consecutive points, big+i^i < ■ ■ ■ < 
big+i.Nig+i j-i, distinct from the singular point 'iTig+i{xo). Actually we have: 

Pig 

1=1 

where mig^ij is the coefficient of the homology matrix Mig. Again, for each j = 
1, . . . we consider the segment: 

(7r/(,+i(a;o), &/o+i,i, . . . , 6;o+i,jv,g+i,j-i, 7r;o+i(xo)). 

we choose the position of the points . . . , 5/(,+i,Ar,^^j^ which minimizes 

the potential energy. Having done it for all loops, we denote Big+i the collec- 
tion of these marked points (all the big+i^^s and 7r;o+i(xo)) on Big+i and con- 
sider the subset of the real line (tt/^+i o I)~^{I3ig^i). It is a discrete subset that 
we can ordered as a bi-infinite increasing sequence (0/q+i,ti)„- This subset of M 
contains the subset {nig+i o I)^^ {nig{xo)) which is a quasicrystal contained in 
the quasicrystal {nig o T)^^{nig{xoj). The configuration (6'/q+i^„)^ is made with 
a concatenation of minimal segments whose extremities are consecutive points in 
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Figure 8. The configurations {Oig^n) and (0io+i,„) 



(tt/q+i o/)~-^(7ri(,_|_i(a;o)) and there are exactly pj^+i equivalence classes of segments, 
each of them corresponding to a minimal segment starting at the beginning and 
ending at the end of a loop in Bi^^+i (See Figure ISJ. 

Step 3: We iterate this procedure to get a configuration {Oig+m.ii)^ foi' each to > 0. 

Lemma 4.3. For each m > 0, the configuration (0;o+m,n)„ has rotation number 
Po- 

Proof. As a preliminary remark, observe that by construction: 

• For any j in {1, . . . ,p{lo)}^ each time the projection of the configuration 
(^io,n)„ crosses the loop jigj of Bi„, it spends an amount of time Ni„,j in 
this loop. 

• Similarly, for any m > and for any k in {1, . . . ,p{lo + fn)}, each time 
the projection of the configuration {Oig+m.n)^ crosses the loop jio+m.k of 
Bi„+m, it spends an amount of time Nig+„i.k in loop. 

• Remark also that for any to > and any k in {1, . . . ,p{lo + rn)}, each 
time the projection of the configuration {Oig^n)^ crosses the loop 7;o+m,fc of 
Big+m, it spends an amount of time in loop which is precisely Ni^+^.k- 

Using the same estimate as for the proof of the existence of a rotation number 
for a minimal configuration, we get that the configuration {9i„^n)n has a rotation 
number and that this rotation number is the limit when to — > +oo of the sequence 
(P™)m>0' where: 

Vm > 0. 



P^0 + "^ 



Claim: The sequence {Pm)m>o constant. 

Proof of the claim: Using the relation (ik-Tkr) we get, for each to > 0: 
1 1 



1 
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Thanks to Lemma 

Plg + m + l 

Thus: 

Prn+l = Prn, Vm > 0. 

This proves the claim and shows that the rotation number of the configuration 
(6'/o.n)„ is equal to po- 

To conclude the proof of the lemma, we remark that a same computation yields 
that, for each p > 0, the configuration (^/g+p n)^ has a rotation number and that 
this rotation number is the limit when m — > +00 of the sequence (Pp,m)„>Q, where: 

Pp.m = iT— Vm > 0. 

J = l 

As shown previously, the sequence iPp,m)m>o constant and Pp.o = Pp = Po- ^ 
Step 4: 

Lemma 4.4. There exists M > such that: 

< 6'(o+m,n+l - Oi„+.m^n < M Vm > 0, \f n ^ Z. 

Proof. Notice first that because of the very construction of the configurations 
(^;o+™,n)n lemma is true if we consider only a finite subset of these sequences. 
Let us prove this lemma by contradiction. Let us fix mp > and assume that 
the lemma is not true for the set of sequences (0io+m,n)„ with m > mg. Choose 
M(mo) > such that M(mo) is larger than the longest loop of Big+m^^. We know 
that there exists m > mo and n G Z such that: 

Af(mo) < 

Recall that the configuration (0/(,+m,n)„ is a concatenation of minimal segments 
whose extremities descend by projection on the singular point of Big+mo- This 
implies that there exists a minimal segment: 

of the configuration (0/o+m,ji)„ ^^'^ ^ loop 7/p+,n,j in Big-^^mo such that: 

o/(e) n7;o+„oj = 0, and 7;o+maj C ni^+,no o I{[9ig+^,n,0i^+^,n+l]). 

Recall that the image Tig^mg^i{jig^mgj) covers all the loops of B/o+mo-i- Using 
Lemma [4. II we deduce that the projection of the segment O on Big+mo-i stays at 
each passage in a loop of Big+mo-i most 3 times in this loop. It follows that the 
rotation number pzo+m of the configuration (0Zo+m,n)„ satisfies: 

1 

Pla+m ^ 77; 7- 
i>l^\.'^la+7no-l) 

This inequality must be true for all mo > and thus po — +00, a contradiction. □ 
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Let US consider the set M equipped with the product topology. For M > 0, the 
set Sm of non decreasing sequences (^n)„ in such that: 

< C„ - ^„-i < M, Vn e Z, 

is a compact subset of M^. Thus it follows from Lemma [4.41 that the set of all the 
configurations (^io+m.n)„7 for m > and their translated is in a compact subset of 

Step 5: For each to > 0, consider Um G M such that belongs to the center of a min- 
imal segment of {Oig^m.n + Um),^- From lemma ETI the sequence of configurations 
(^io+m,n + Um),j has an accumulation point in M^. We denote this configuration 

(^00,Tl)„- 

Lemma 4.5. The configuration (^?oo,n)„ *s a minimal configuration with rotation 
number pQ. 

Proof. The fact that the configuration (0oo,n)„ is minimal is standard. Consider a 
segment of (0oo,n)„- By construction this segment is a limit of minimal segments 
and it is straightforward to show that this segment is minimal. 
Let us prove now that the configuration (^oo,ri)„ has rotation number pQ. Since the 
configuration is minimal, it has a rotation number which is defined as the limit: 

i=i 

where Noo,i,j is the minimal number of times the configuration (0oo.n)„ spends in 
the loop oiBi. 

We use a similar argument to the one used in the proof of the continuity of the 
rotation number. Fix li > Iq, and choose a loop ji-^j in Bi-^. Consider the first 
time when, starting from on the real line and going in the positive direction, the 
configuration (0oo,n)„ enters in this loop. Let us do the same for the configuration 
{Oi„+m,n + Wm)„- Since a subsequence of configurations (0i(,+m,n + 'i'in)„ converges, 
when TO goes +00 to the configuration {6oo,n)„, it follows that for to big enough 
both projections of the configurations stay the same time in the loop ji-^j for their 
first visit in this loop. It follows from Lemma l4.1l that the minimal number Ni^-^-m.i.j 
of times the projection of the configuration {9ig+m,n)n spends in the loop 7;^^- of 
Bi^ satisfies: 

- Noo,i,,j\ < 2, Vj e {!,..., pGi)}. 



Thus for TO big enough, the rotation number po of the configuration 
satisfies: 

1 1 



< PO < 



Pli — ' " — Pli 



On the other hand 



1 1 



Pll _ roc, _ p,^ 

E Vh,jNoo,h,j + 2t/(C;J Vh,]Noc,h,j 
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This implies: 

1 1 



Poo pa 

Since this last inequality is true for alHi > we get: 

Poo = Po- 

□ 

This ends the proof of Proposition 14.21 □ 

In order to prove Part [iii) of Theorem II .11 we choose a positive real number p 
and consider a sequence of minimal configurations {Om.n)^^ fn > 0, with rotation 
number pm G ■J' such that: 

lim Pm = p- 

m — > + oo 

A discussion completely similar to the one we used in the proof of Lemma l4 . 41 allows 
us to show that there exists M > such that: 

< 9„,,7i+i - 0,nji < M Vm > 0, Vn e Z. 

Consequently, the set of all the configurations {9rn.n)m for to > and their trans- 
lated, is in a compact subset of and thus, as done previously, we can exhibit 
a subsequence of configurations which converges to a minimal configuration {0n)n- 
Thanks to continuity property of the rotation number (Part (ii) of Theorem II. 
we conclude that the rotation number of (0n)„ is p- 

5. Final remarks 

5.1. Dynamical systems. Minimal configurations of the Frenkel-Kontorova model 
obviously satisfy the variational equations: 

U'{e^ - On+i) - U'iOn-i - e„) + V\0n) = 0, Vn S Z. 

By introducing the new variables^: 

p„ = [/'(0„_i - 0„), VneZ, 

we get the dynamical system defined on M x M by: 

Pn+l ^ Pn - V'{6n) 

(★ ★ ★) 

On+1 = On - {U')-\Pn ~ ^'(^n)) 

In the crystal case, V' is a periodic function with period L, the period of the crys- 
tal. It follows that the map defined by (★★★) descends to a map on the open annulus 
R/L.R X R which is an orientation preserving diffeomorphism which preserves the 
standard area form. Area preserving maps of the annulus have been widely studied 
and Aubry-Mather theory which makes a bridge between the Frenkel-Kontorova 
model and dynamical systems, has been a powerful tool for both sides. 

In the quasicrystal case, the dynamical system extends to an area preserving 
" diffeomomorphism" ® on the solenoidal annulus Vl[QC) x M. The study of such 
maps will be the subject of a forthcoming paper. 



^Recall that U' is an increasing homeomorphism of the real hne. 
^We mean difEeomomorphism in the leaf direction. 
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5.2. Quasicrystals in R'^, d > 1. As we already noticed, the construction of 
the hull of a quasicrystal and its interpretation as an inverse limit of branched 
manifolds can be done for quasicrystals in any dimension (see |EG) . |BBGj . 
On the other hand, in a recent work KLR^ , H. Koch, R de la Llave and C. Radin 
developed a generalization of Aubry-Mather theory for functions on lattices in R'*. 
Both arguments make tempting to develop in a same way, a Aubry-Mather theory 
for quasicrystals in M'^, d> 1. 
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